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The dynamics of neutrino mixing and oscillations are studied directly in finite real time 
in a model that effectively describes charged current weak interactions. Finite time cor- 
rections to the S-matrix result for the appearance and disappearance probabilities are 
obtained. It is observed that these effects may be of the same order of the S-matrix result 
in long-baseline appearance experiments. We argue that fundamentally, the S-matrix is 
ill-suited to describe long-baseline events due to the fact that the neutrino is produced in 
an entangled state with the charged lepton, which can be disentangled by the measure- 
ment of the charged lepton near the production site. The appearance and disappearance 
far-detection process is described from the time evolution of this disentangled "collapsed" 
state, allowing us to establish the conditions under which factorization of detection rates 
emerges in long-baseline experiments. We also study the time evolution of the reduced 
density matrix and show explicitly how oscillations are manifest in the off-diagonal terms, 
i.e., coherences, as a result of a finite time analysis. Lastly, we study a model for the 
"GSI anomaly" obtaining the time evolution of the population of parent and daughter 
particles directly in real time. We confirm that the decay rate of parent and growth rate 
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of daughters do not feature oscillatory behavior from interference of mass eigenstatcs. 
Keywords: neutrinos oscillation; quantum field theory; entanglement. 
PACS numbers: 14.60.Pq;13.15.+g;12.15.Ff 

1. Introduction 

Neutrino masses, mixing and oscillati ons are the clearest evidence yet of 
physics beyond the standa rd mod el I ^P I They provide an explanation of 
the solar neutrino problem I 4 | 5 | 6 | 7 | an( j have oth er important phenomenologi- 
cal l l | 2 | 3 | 8 | 9 | 10 | ll | 12 | 13 | 14 | 15 ] astrophysical I 5 | 16 | 17 | anc j cosmo logical I 18 | 19 | conse- 
quences. Moreover, another fascinating aspect of mixing and oscillations of ultrarel- 
ativistic neutrinos is that they provide a remarkable manifestation of macroscopic 
quantum coherence over unprecedented long distance and time scales. Whereas in 
condensed matter systems macroscopic quantum coherence is maintained over meso- 
scopic scales, long baseline disappearance and appearance experiments probe the 
coherence of the neutrino states over scales of hundreds of kilometers. 

In its simplest (and perhaps overly naive) inception, mixing and oscillations 
of massive neutrinos have been explained by an analogy with a two level system 
undergoing Rabi-like oscillations between them. Oscillations emerge from the in- 
terference of the quantum m echanical states associated with the mass eigenstates 
(see for example | 1 | 2 | 3 | 9 | 10 J8| an( j re f erenC es therein) . As simple and compelling as 
this interpretation is, deeper investigations of t his basic paradigm have raised a 
number of important and fundame ntal quest ions j^j. O ne of these involves the en- 
ergy and mo mentum unc ertainties | 22 | 23 | 24 | 25 | 26 | 27 j28| a ^ pj c ^ n at is still receiv- 
ing attention | 29 | 30 | 31 | 32 | ^ or a recen t rev iew see ^^j). The recognition that exact 
energy and moment um conservation prevent oscillations between neutrinos of dif- 

lOOIOQI 

ferent masses 1 1 1 has led to the consideration of oscillation experiments in terms 
of wave-packets | 22 | 23 | 24 |25| j nc luding the quantum mechanical aspects of produc- 
tion and d etection which have been incorporated in the quantum field t heoretical 
framework I^4i:i^i:ffil37l38l39l4()l41l42l43l44l45l46l47l48l49l5()l51l52l53l54l33l 

Previous quantum field theory treatments of mixing and oscillations are 

S-matrix theoretic in nature, making use of in-out wave-packets spatially 

localized at the source, or the "near" detecto jf], and the far detector 
I:ffil37l38l39l4()l41l42l43l44l45l46l47l48l49l5()l51l52l54l33l Howeverj in all of these 

treatments the S-matrix calculation takes the interaction time to infinity, even when 

a We are using the term "near" detector to describe measurement apparatus at or nearby the 
production site, which is used to detect the charged lepton that is produced with the primary 
neutrino. This is not to be confused wit h t he term near detector that is widely used in the 
experimental literatures (see for example ESI) . The latter is used to detect neutrino at a short 
baseline. 



January 26, 2012 1:10 WSPC/INSTRUCTION FILE nuoscillation-IJMPA- 
revised 



Neutrino Oscillations, Entanglement and Coherence 3 

in some treatment s the initial an d final wave packets are denned at finite source 
and detector times HS )50 | 51 | 52 |54| j n wave-packet-S-matrix approach advocated 
in many of the references above, the wave packets spatially localized at the source 
and far detectors respectively are spatially separated by the finite baseline, and also 
localized in time, however, the intermediate steps invoke the S-matrix approach 
wherein in and out states are prepared at fcj = — oo,iy=oo . 

This incongruity between keeping a finite baseline, with wave packets defined at 
some initial time at the source and final time at the far detector, and taking the 
time to infinity in the S-matrix element is usually justified with the statement that 
time is not measured in appearance or disappearance experiments. 

While this is true, it is a practical, but not a fundamental reason for taking 
the infinite time limit. One can consider a gedanken experiment in which clocks at 
the source (near detector) and far detector are synchronized via global positioning 
satellites and register the detection of the charged leptons at the source and far de- 
tectors at the time ts and tp, respectively. Obviously, the time difference registered 
by these clocks trj — t$ ~ L, with L being the baseline. This has uncertainties of 
the order of the size of the source and detectors as the interaction vertices from 
which the charged leptons emerge are localized within these regions (although cur- 
rent resolution of the interaction vertices is much more accurate than this scale) . 

Furthermore, taking the infinite time limit as in the usual S-matrix calculation 
enforces total energy conservation via an overall delta function in the transition am- 
plitude. The transition probability treats the square of this delta function as overall 
energy conservation multiplied by the total time of the interaction, from which an 
interaction rate is extracted by dividing by this time (in the long time limit). How- 
ever, oscillations arising from quantum mechanical interference have nothing to do 
with a transition rate and in principle do not feature a secular evolution in time. 

While there are various quantum field theory calculations of appearance and dis- 
appearance processes in the literature (see references above), all of them invoke the 
S-matrix approach and take the infinite time limit from the outset. It is fundamen- 
tally important to understand any possible caveat in S-matrix approach by focusing 
on the fact that only a finite time interval elapses between the production and the 
detection of neutrinos. In this article, we aim to demonstrate that conventional S- 
matrix calculation where in and out states are eigenstates of energy-momentum is 
ill-suited in principle to describe neutrino oscillations by studying the quantum field 
theory of neutrino mixing and oscillations in real time in the plane wave limit. We 
keep a finite time interval between the initial state "prepared" at the source and the 
final state measured at the (far) detector. After obtaining the transition amplitudes 
and probabilities, we point out the problems in the standard S-matrix calculation. 
Then, combined wit h ou r recent work on the disentanglement issue associated with 
neutrino oscillations we seek a complete space-time description of the dynam- 
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ics of mixing and oscillations that describes long-baseline experiments in the plane 
wave limit. 

We obtain the transition matrix element in which neutrinos are produced and 
detected via charged "leptons" , directly in real time. Comparing with the S-matrix 
results in the plane wave basis, we discuss how wave packet localization restricts 
the contributions from different channels and obtain the non-conventional transi- 
tion probabilities at finite but long time. For the appearance probability, we find that 
finite time corrections are comparable to or larger than the S-matrix contribution. 

Then we argue that this formulation is incorrect because long-baseline exper- 
iments involve detection measurements at two different times, corresponding to 
detection or "collapse" of the charged leptons at the near and the far detectors. 
The neutrino is produced in an entangled state at the interaction vertex at the 
source, but the measurement of the charged lepton at the near detector disentan- 
gles the neutrino. It is this "collapsed" disentangled state that evolves in time and 
leads to the production of the charged lepton which is measured in the far detector 
| 46 | 32 | 56 [57| carr y ou ^ a systematical study containing the dynamics of entangle- 
ment and disentanglement, namely that neutrinos produced by the decay of a parent 
particle are entangled with the daughter charged le pton a nd then disentangled close 
to the near detector. Different from the treatment | 46 | 32 | mvo king the infinite time 
limit, we analyze the time evolution of the reduced density matrix for the neutrino 
and show that keeping finite time allows us to understand the time evolution of the 
population and coherences. We establish that indeed coherence, as determined by 
the off-diagonal dens ity matrix elements in the mass basis is maintained up to the 
oscillation time scale 

We then obtain the near and far detector event rates from the detailed evolution 
of the entangled state, disentanglement at the near detector and further evolution 
to the far detector. We discuss under which circumstances the factorization of the 
processes is valid. We also discuss in detail under what circumstances the usual 
quantum mechanical description is valid. The real-time analysis of disentanglement 
and coherence shows that this is the case provided the neutrino state is disentan- 
gled on time scales much shorter than the oscillation scale. This analysis also makes 
manifest the compatibility with energy conservation. 

Another clear evidence of the fact that the theory of neutrino oscillations is 
far from being understood at a fundamental level is found in recent controver- 
sies re garding the rec oilless emission and detection of neutrinos (Mossbauer neu- 
trinos )TSHl59l60l6ll62^Y(^ found in the controversial interpretation 
of the GSI anomaly I 63 | 64 | 31 | 65 | 66 | 67 86 name iy ; periodic modulations in the K- 
electron capture and /3 + decay rate I ' ' in terms of quantum beats resulting from 
the mixing of neutrino mass eigenstates and their interference in the final state 
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| 70 | 71 | 72 | 73 | 74 | 86 ] j n ^ g ar -ticle, we review our real time analysis of the decay of 
a parent particle into a daughter particle and a "flavor" neutrino ^3 as an appli- 
cation of the real time method developed for neutrino oscillations. We show that 
both the decay rate of the parent, and growth rate of the daughter populations do 
not feature oscillations as a consequence of mixing and establish unambiguously the 
reasons. We therefore conclude that the "GSI-anomaly" cannot be explained by the 
interference of mass eigenstates in the final state. 



2. A Model of "Neutrino" Oscillations 

The goal of this work is to study the dynamics of mixing and oscillation of neu- 
trinos in quantum field theory but in a finite time interval. In order to exhibit the 
main results in a clear and simpler manner, we introduce a bosonic model that de- 
scribes mixing, oscillations and charged current weak interactions reliably without 
the complications associated with fermionic and gauge fields. We can do so because 
the technical complications associated with spinors and gauge fields are irrelevant 
to the physics of mixing and oscillations, as is obviously manifest in meson mixing. 
Our model is defined by the following Lagrangian density 

C = Co[W, I a ] + C [v a ] +£int[W, l a ,Va] 5 Oi = e, /J (1) 

with 

C [v] = ~ [d^ T d^ - f T Mf] , (2) 
where ^ is a flavor doublet representing the neutrinos 

*=(:)• < 3 » 



and M is the mass matrix 



(4) 



The interaction Lagrangian is similar to the charged current interaction of the stan- 
dard model, namely 

C- mt (x,t) = gW(x,t) l e (x,t) v e (x,t) + ln(x,t) v^{x,t) , (5) 

where g is the coupling constant. W(x) represents the vector boson, or alternatively 
the pion field, and l a , a — e, \i the two charged leptons. From the interaction 
Lagrangian (JSj) , it is clear that W, l a and v a are all hermitian. The mass matrix is 
diagonalized by a unitary transformation 

u-\e)Mu(9) = ( m n l0 ) ; u{6) = ( cos \ sin ^) . (6) 

\ 77i 2 / Y — sin cos J 
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In terms of the doublet of mass eigenstates, the flavor doublet can be expressed as 



This bosonic model clearly describes charged current weak interactions reliably as 
it includes all the relevant aspects of mixing and oscillations. 

We consider "neutrino" oscillation experiments following the interaction pro- 
cesses illustrated in Fig. [U 

W i + v f W + l P 'P^a appearance 

1 W + l a disappearance ' ' 





Fig. 1. Typical experiment in which the charged leptons are measured at a near and far detector 
and the neutrino is an intermediate state. 



In a quantum field theory calculation of the transition amplitudes, neutrinos 
propagating betw een the production and detection regions are considered as in ter- 
mediate particles I24I34I35I36I37I38I39I40I41I42I43I44I45I46I50I51I52I48I47I54I33I and 

are described by a (free field) propagator. 

In order to extract the main physical aspects and new results in the clearest 
manner, we study the time evolution in terms of plane waves, postponing the nec- 
essarily more technical discussion in terms of wave-packets to a follow up article. 
However, to set the stage to the next step that includes such formulation, we revisit 
some aspects related to wave packet localization that are relevant for the present 
discussion. 

For a wave packet description we consider the initial and final W, states to be 
described by wave packets. Specifically, we consider the following type of spatially 
localized states 

f(x,X;p) = N e- {JL ^£ e - ** (9) 

where N is a normalization factor, X$,d is the center of the wave packet for source 
and detection, respectively, a is the localization length of the wave packet, which 
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might be a nuclear scale or a macroscopic scale of the order of the size of the 
detector and production regions. We will only require that a <C \Xd — Xs\- The 
spatial Fourier transform of these wave packets is 



d 3 xe* p - s f(x, X-p)=Af e-W-ti 1 ^ e H?H0-* . = 



'2tz 



(10) 

The normalization factor in Fourier space, Af, has been chosen so that in the 
macroscopic limit a — > oo, the wave packet has a definite momentum f(p,k;X) — > 
fc). The initial and final quantum states are localized at the source (Xs) or 
detection (Xp) sites and are described by wave packets that feature mean momenta 
K S ;K D ;P S , P D respectively 



\W(S)) = 


= \W(K s ;Xs)) = 


f dk s f(Ks,k s ;X s ) \W(k s )), 




\W(D)) = 


\W(K D ;X D )) = 


1 dk D f(K D ,k D ;X D ) \W(k D )), 




MS)) 


= \T a (Ps\Xs)) = 


f dps f(P s ,Ps;Xs) \la(ps)), 




\T a (D)) = 


ee \l a (P D ;X D )) = 


1 dp D f(P D ,p D ;X D ) \L(pd)), 


(11) 



where the quantum states \W(ks)) , \W(po)) , \l a (ks)) > |£a(Pz>)) on the right hand 
side of Eq. (jllj) are plane wave single particle states, which are the eigenstates of 
the non-interacting Hamiltonian. In term of these localized states, the appearance 
and disappearance transition amplitudes for wave packets prepared and detected 
at the source (near detector) at an initial time fcj and detected at a far detector at 
a final time tf are given by 

A a ^ = (W(D)MD);T a (S)\e- ift ^-^\W(S)), (12) 
A Q ^ a = (W{D)MD)-MS)\e- i6 ^-^\W(S)), (13) 

where D and S label the states localized in the detection and source regions, re- 
spectively. The respective probabilities are 

V a ^ fj = \A a ^p\ 2 and V a ^ a = \A a ^ a \ 2 . (14) 

The transition amplitudes between the initial and final localized states are given 

by 

(W(D)MD)-MS)\e'^^\W(S)) 
= J dk s dps dp D dk D e-^f A a ^ e iBit * f*(K D ,k D ;X D ) 
f*(Ps,Ps; X s ) f*(P D ,PD;X D ) f(K s , k s ; Xs), (15) 
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where 

E l = El ; E f = El+El D +El s , (16) 
and ti and tf label the initial and final times. 

Aa^p = (W(k D ),l {p D );l a (p s )\U(t f ,ti)\W(ks)) (17) 

is the usual transition matrix element in terms of plane waves eigenstates of the non- 
interacting Hamiltonian in the interaction picture. To obtain the above expressions 
we have passed to the interaction picture by writing 

e -iA(t f - ti ) = e -iH t f u(t hU ) e iH t t ; ( lg ) 

with the time evolution operator in the interaction picture 

U(f f ,ti) = e lHot f e -&(*f-*i) e~ lHoU = T e l/ ''/ d ' 3xdt c "^ s ^ , (19) 

and we have also used the fact that the initial and final quantum states in (|17[) are 
free single particle plane wave eigenstates of the non-interacting Hamiltonian Hq 
with Ef = y/tf+Mjy, El = y/p^+mf. 

These are the main ingredients in the appearance and disappearance transition 
probabilities in terms of the localized wave packets. In the limit when the localization 
length goes to infinity a — > oo with the normalization (jlOp . the functions / become 
delta functions in momentum and the localized states become simple plane wave 
states. In what follows we will obtain the transition matrix element A a ^/3 (|17p up 
to second order in the interaction. 



3. Appearance and disappearance amplitudes and probabilities 

3.1. Localization and suppression of crossed channel contributions 

Up to second order in the coupling g the matrix element A a ^p features the contri- 
butions depicted in figs. (|2^,b) respectively. If the lepton at the production vertex 
were in-coming, Fig. (a) would be the equivalent of an s-channel and Fig. (b) of 
a t-channel contribution. We refer to these as "s" and "t" channel contributions 
respectively for the remainder of the discussion. 

For the "s" and "t" -channels we find respectively 

A^lp = -g 2 n J*' dhdt 2 J d 3 Xl d 3 x 2 e i ^ tl -^- ari ) e i(^r fl ti-fe-irO e <(4 s *»-ft-&) 

e -i{EY s t 2 -k s .S 2 ) (olT^OnHte)) |0), (20) 

and 

A { ^ = -9 2 n J*' dhdt 2 J d 3 Xl d 3 x 2 e^^-^^Ogi^fa-fo-aOe^ta-fo-fe) 

e-'Ws**-**- 3 *) (0\T(^( Xl )u a (x 2 ))\0), (21) 
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(a) 

Fig. 2. Feynman diagrams that contribute to _4 Q 
to the "s-channel" and (b) to the "t-channel." 




up to second order. Diagram (a) corresponds 



where for disappearance a — /3, for appearance a^(3, and 

1 



n 



16 V 4 EW El El 

k s k D PS PD 



(22) 



The transition matrix elements A a ->p between the initial and final localized 
states are obtained from the A a ^>/3 matrix elements above by convolution with the 
wave packets describing these localized states (see Eq. [15)1 . In order to understand 
the features associated with the localization, we consider wave packets with a macro- 
scopic localization scale a <C \Xd — Xs\ of the order of the size of the detectors, 
in which case the wave packets are nearly plane waves and we can approximate 
ps,D — Ps,D, ks,D — Ks,d in the arguments. For the s-channel we obtain 

(X!-X S ) 2 



^ 



whereas for the t-channel contribution we obtain 



7i — 



2.2 



— ^ — 



(g2-£pr 

^2 



(23) 



(24) 



Therefore, while for the s-channel the x\,x% integrals factor out, the t-channel con- 
tribution obviously vanishes for \Xd — Xs\ S> er. In other words, a localization of 
the initial and final states even over macroscopic scales, as long as a <C \Xd — Xg\, 
ensures that the "t" -channel contribution depicted in the second Feynman diagram 
Fig. [2] vanishes. 

3.2. Amplitudes and probabilities for plane waves 

In what follows we consider only the "s" -channel contribution and focus solely on 
studying the transition amplitudes for plane waves, postponing to forthcoming work 
a full space-time description of the production-detection process. 



The plane wave transition amplitudes A a ~^p are given by 



A a ^f3 = —g II J dt\dt2 I d x\d X2 e 



i{E% ti—kn-Si) e i(E l D ti-fin-xi) e i(E l s t 2 -p s -S 2 ) 

e - l{ E^t 2 -ks-s 2 ) ( \T( V p( Xl )v a (x2))\0)(25) 



January 26, 2012 1:10 WSPC/INSTRUCTION FILE nuoscillation-IJMPA- 
revised 



10 Jun Wu, Jimmy A. Hutasoit, Daniel Boyanovsky, Richard Holman 

Here 

(0\T(ve( Xl )v a (x2))\0) = UpjUjaiOlTfoixMxtWO) , (26) 

J'=l,2 

where U a b are the matrix elements of the matrix U given in Eq. ([6]), and 

are the propagators for the mass eigenstates, with 

0| = p 2 + m 2 . (28) 

Here, 25]j (E l D ) and (p*o) label the energy and momentum of the vector boson 
(charged lepton) at the detection region, while Eg {E l s ) and k s (ps) label those of 
the production region. 

The spatial integrals over X\ and X2 yield the usual spatial momentum conserva- 
tion (27r) 3 <5(p — k D —po) and (27r) 3 (5(fcg — p s ~ p), respectively, leading to an overall 
momentum conservation in the amplitudes A a ->/3 for plane waves. In the transition 
amplitude for wave packets total momentum conservation is smeared by the wave 
packet extension. 

The integrals over the finite time intervals 

/ = / / dhe^-"^ f 1 dt 2 e-^ Es ~^ ■ Es — Eg — E l ffs ; E D = Eg + El D 

Jti Jti S D • - 

(29) 

require a careful treatment. In S-matrix theory, ti — > —oo ; tf — > oo and the integrals 
require an adiabatic switching-on convergence factor. We want to analyze this long- 
time limit to establish contact with the S-matrix results. It proves convenient to 
write 

T t T t 

*/=2 + 2' k= 2~2 (30) 
and to introduce an adiabatic convergence factor in the integrals 

t t 
1= lim e lAjE * /* a e <(-Bx>-«)*i e -*l*il ^ [' e -i(E s -u)t a e -s\u\ ^ ( 31 ) 

where 

AE = E D -E s = E f -Ei, (32) 

and the initial and final energies, Ei and Ef, are defined in Eq. (fT6|) . We find 

t 

J 2 e i{E-u)H e -S\H\ dh =27r S ( E _ U -) + 2 S [(E-u);t , (33) 
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where the function (distribution) 



S 



iH0+ 



(E-w);t = - lim / cos [(E - u)ti] e~ sltl1 dt 1 



has the following behavior 
(E-u);t 



lim S 

t— >oo 



; lim S 



{E-u);t =-tt8{E-w). 



(34) 



(35) 



For a finite time interval and taking first the limit S — > + , it follows from (|3"31 that 



7T 5(E-lu) +S (E-uj);t 



(5^0+ 



We also gather the following useful results 



lim 



lim 



<5^o+ i(uj — E + i5) 5->o+ 



sin 




-0,)§ 




(E- 






-E4 





(36) 



lim 



= — lim 

(5^o+ i(uj — E — id) 5^o+ 



(37) 



which both go to zero at the infinite time limit t — > oo. 



We note that the prefactor e lAE ¥ in (|3"Tj) combines with the factor e~ l ^ E f t t~ Eiti ^ 
in (fl~5f to yield e~ lE ( t f~ ti '> where -E = (Ef+Ei)/2, manifestly displaying time trans- 
lational invariance. 



The plane wave transition amplitudes are then given by 

Aa^p{t) = -g 2 n (2ir) 3 S(k s - p s - k D - p D ) Up 3 l[fi i; t] U ja , (38) 
where X \flj ; t\ is given by the dispersive integral 

{tt«5(w - E s ) +S[(u - E 8 )]t] }{tt«5(w - £ D ) + «s[(cj - E D );tj } 



lu 2 - n 2 + it 



(39) 



and the momentum argument of the frequencies fl,- is p = ko + Pd = ks — ps- 

Before we carry out the u) integral, we note that taking the t — > oo limit at this 
stage yields energy-momentum conservation at each vertex. This yields 

2n5(Ef -E{) 



El - fi? + ie 



(40) 



J s --j 

leading to the usual S-matrix result for the transition amplitude for plane wave 
initial and final states 

' : U jai (41) 



A^-J^ = -g 2 U(2n) 4 S(E f -E i ) 6(k s -p s -k D ~p D ) U Pj 



e 2 - n 2 + ,< 
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clearly demonstrating that no oscillations occur for plane waves in the infinite time 
limit. 

For finite time interval, the integration over lo (j39]) is tedious but straightforward. 
The result is 



„— if2.it 



e% - nj + it 



where to simplify notation we have introduced 

D(E f - Ei]t) = 2-K§(Ef - Ei) + 2S {E f - Ei);t 



(42) 



(43) 



(E f -Ei+2i8)± 



2Qj 1 (Ef - Ei + 2iS) 

e -i{E f -Ei-2i5)i 



E f - Ei 



(E f - Ei - 2iS) 



(Clj + E s )(Clj +E D ) 

Ef — Ej 

(Slj -EsWj -E D ) 



(44) 



and 



, , e -i(E n +Es-2iS)± 

2nA {E D + E s - 2iS) 



Ed + E s 



e i(E D +E s +2i8)± 

' (E D +E s + 2iS) 



(Qj +E s )(n j +E D ) 
Ed + Es 



(Slj - E s m - E D ) 



(45) 



In the t — y oo limit, the first term in (|42[) gives the S-matrix result. The functions 
Ji and T only depend on the frequencies flj in the denominators and vanish in the 
t — y oo limit, as can be seen from Eq. (f3"T|) . It is straightforward to show that there 
are no poles at Qj — ±Ed,s because the residues vanish, so that both % and T 
must be understood in terms of their principal part. 

The disappearance transition amplitude for plane waves are then given by 



-g 2 n (2nf5(k s -ps-k D - p D ) 



cos 2 9 T[Qi;t] + sin 2 9 l[fl 2 ; t] 



= -9 2 n (2n) 3 S(k s -ps - k D -p D ) 



sin 2 9l[n 1 ;t]+ cos 2 9 X[Q 2 ;i] 



while the appearance amplitude is given by 
A e ^^ = g 2 II (2ir) 3 5(k s - p s - k D - p D ) cosd sin( 



x[fii;t] -x[n 3 ;t] 



(46) 



(47) 



(48) 



Taking the limit t — > oo in these amplitudes, one recovers the S-matrix result 
obtained by replacing 2Tfij; oo] by (|40l) . in which case there is no time dependence 



January 26, 2012 1:10 WSPC/INSTRUCTION FILE nuoscillation-IJMPA- 
rcvised 



Neutrino Oscillations, Entanglement and Coherence 13 

and no oscillations from interference terms. Instead of taking this limit, we consider 
a time interval large compared to microscopic times but of the order of the oscilla- 
tion time corresponding to the experimental situation in which the baseline is just 
long enough for a few oscillations. For finite time, we can set 5 — > + in the above 
expressions, since we have explicitly separated the delta functions and identified 
which terms vanish in the formal limit t — > oo. 

We would like to emphasize that there is an important distinction between what 
we do here and the usual S-matrix approach. In the S-matrix approach, the tran- 
sition amplitude is obtained in the t — > oo limit and after taking this limit, one 
obtains the transition probability. In contrast, we obtain the transition amplitude 
at finite time t and obtain the probability. 

In order to extract the most relevant contributions at long time, we consider the 
ultrarclativistic limit and write 



q 1 = q-a, n 2 = n + A, 



(49) 



where 



n = 



o m? + m% 
P + — - 



a c 2 2 2 

A = —=-, dm =m 9 — m-,, 

4fl 2 1, 



(50) 



taking A <C ^ as is the case for ultrarclativistic neutrinos. We can then write 



W[fij-;t] =u\n-t\ +(-iy 



dn[fl;t] 
dTt 




(51) 



We consider the realistic situation in which Sm 2 /fl <C 1 and keep the terms of 
0(Sm 2 /Q ) only in the S-matrix term and in the exponentials e* 012 *, neglecting 
the small corrections both in W and T. The reason for keeping the 0(Sm 2 /Q 2 ) 
correction in the S-matrix contribution will become clear below when we discuss 
the appearance probability. 



We can further define 

F[U;t] =e- tTit F\tt-i\, (52) 
thus simplifying the expressions for the disappearance and appearance amplitudes 

cos 2 f) sin 2 f) 

A e ^ e = -g A Il{2TTf5{k s -p s -k D -p D ){D{E f -E l -t) 



+ 



E 2 s -n\ + ie E 2 -n 2 + it 



+ n\Q;t] +F[Q;t] 



cos de *» + sin e 



O 



Sm 2 



(53) 
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A e ^ = g 2 n (2tt) 3 S(k s - p s -k D -p D ) cosB sin6 { D(E f - E % :t) 



■2iF\n-,t] sin 



Sm 2 ' 
4ft ■ 



O 



ft 2 



El 



ft 2 



El 



(54) 



The amplitude A^^ can be obtained from (|53| by the replacement sin 2 9 f-> cos 2 9. 

The expressions above clearly exhibit how and where oscillatory interference 
terms arise in the probabilities. The usual S-matrix result is obtained in the t — > oo 
limit where D(Ef — Ef, t) -> 2w5(Ef — Ei) and F\Q, t], Hp, t] -> 0. It is clear that 
in this limit, the oscillatory behavior is suppressed and no interference terms can 
possibly survive in the transition probabilities. 

In obtaining the transition probabilities, we recognize two types of oscillatory 
terms: terms that feature exponentials of the form e ± l ( E o±Es)t anc j ^ nose f the 
form e ±iSm */ 4 ". The former are fast oscillatory terms on microscopic time scales, 
whereas the latter are slow phases on these time scales and only manifest themselves 
on much longer time scales, of order of the baseline t ~ L. The fast oscillatory 
terms average out on the (much) longer time scale. After obtaining the transition 
probabilities at a finite time t, we discard terms that feature the fast oscillations 
that average out in the long time li mit, akin to what happens in the "rotating wave 
approximation" in quantum optics and finally take the 5 — > limit, obtaining 
the plane wave transition probabilities 



g 4 n 2 V 6(ks - ps - k D - p D )<2TTt S(E f - Ei) 



cos' 



sm 



Ei-n 2 



+ 



+ 



4ft 2 



1 



(n + E s ) 2 (n + E D ) 2 {n-E s ) 2 {n-E D ) 2 
i i 



4 rr (n + E s ) 2 (n + E D ) 2 (n-E s ) 2 (n-E D ) 2 

g 4 LT 2 V (2Trf5{k s -p s -k D ~ p D ) sin 2 26 

x \ t 6(Er - Ei) 3_ L 

I 4 V / ' E 2 s -n 2 + ie El~n 2 + ie 

1 1 



4ft 



(n + E s ) 2 {n + E D ) 2 (n-E s ) 2 (n~E D ) 2 



E 2 -n 2 



J S "1 T ^ J^g 

1 - sin 2 29 sin 2 



Sm 2 ' 
Iff 



4ft 



(55) 



(56) 



Here, we have neglected contributions of 0(5m 2 /ft ). As before, the disappearance 
probability is obtained from (|55ll by the substitution cos 2 9 O sin 2 9. 
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The first line in the above expressions is the S-matrix result, where we have used 
(2ir8{Ef — Ei)) 2 — » (2n) t S(Ef — Ei) as usual. The second line in (|55|) arises from 
|%(f2;£)| . This term is a direct result of calculating the probability at finite time 
before taking the t — > oo limit. The two procedures do not commute. On one hand, 
taking the t — > oo limit first results in the vanishing of % by the averaging of the 
oscillatory terms (Riemann-Lebesgue lemma). On the other hand, obtaining the 
probability first and taking the long time limit afterward yields the contribution 
from the modulus squared of each oscillatory term, leading to the second line in 
IB- 



Writing as in Eqs. (|49I50[) . it follows that the S-matrix contribution to the 
appearance probability is 



V^T « m - Ek) (57) 



An 



Upon integrating on the final density of states for times such that there are few 
oscillations, namely Sm 2 t/4fl ~ 1, the S-matrix contribution is much smaller than 
the oscillatory terms in the second line of ((56)) because 5m 2 <C ,Eg D . On the 
other hand, the S-matrix contribution to the disappearance probability is 

V s e -T tx <xt6(E,- E^ l _ 2 , (58) 
(E 2 - n y 

which upon integration over the density of states will dominate over the oscillatory 
terms in these probabilities for t ~ t osc ~ Afl/Sm 2 . Thus we see that appearance 
and disappearance probabilities are fundamentally different. For time scales during 
which oscillation phenomena can be observed, the S-matrix contribution to the dis- 
appearance probability dominates while this may not be the case for the appearance 
probability, as the oscillatory terms emerging at finite time may be comparable to 
or larger than the S-matrix contributions. 

Lastly, we would also like to note that the oscillatory terms feature an energy 
dependence very different from that of the S-matrix contribution. 

3.3. Which time scale? 

An important question emerges from the analysis above. What is the value of the 
time t in these probabilities given in Eqs. i55\) and [5b]) ? The linear time dependence, 
a consequenc e of total energy conservation, also emerges in many other a nalyses with 
wave padxtG 123 ™ 71381 ^' 401 ^^ from which 

the S-matrix evaluation of the transition probabilities always yields a result pro- 
portional to T/L 2 where T is the "total reaction time" arising from squaring the 
delta function associated with overall energy conservation and L is the baseline. The 
factor 1 / L 2 has a clear physical meaning corresponding to the ratio of the neutrino 
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flux received at the far detector and that produced at the near detector. In the plane 
wave limit, the factor 1/L 2 disappears because we lose the localization property of 
wave packets; however, the time t in the probabilities (|55[) and (|56[) plays the same 
role as the linear time dependence T in standard S-matrix calculations. This is clear 
shown in Eqs. (|57[) and (|58p denoting the S-matrix contributions to the probabilities 
given in t|55|) and (|56[) . However, the finite time corrections to the appearance and 
disappearance probabilities in (|55[) and (1561) do not feature the linear time depen- 
dence. This observation exactly motivates us to address the question what the time 
scale t really means. 

In a typical S-matrix calculation one divides by T, the "total reaction time," 
to obtain a transition rate. However, this interpretation needs revision in the case 
of long-baseline experiments. In these experiments, a charged lepton is measured 
at the near detector (source) whereas another charged lepton is measured at the 
far detector. The intermediate neutrino state propagates between these as a wave 
packet. Therefore, there are two time scales in this case: the time at which the 
near-detector measurement of the charged lepton occurs and the time at which the 
detection of the charged lepton at the far detector occurs. The question of what 
precisely is the time t in the S-matrix calculation is in principle independent from 
the wave packet treatment and is an inherent question to the S-matrix formulation 
of the production, detection and propagation in long-baseline experiments. 

There is another problem with the analysis that we have done so far. In the 
absence of oscillations, dividing by time or taking the time derivative yields identi- 
cal results. However, in the presence of oscillatory contributions, the transition rate 
must be obtained by taking the time derivative. The time derivative of the oscilla- 
tory terms featured in the transition probabilities (|55H56p are of order 5m 2 /CI. For 
the disappearance probabilities, this derivative term is subleading with respect to 
the S-matrix contribution, but it may be of the same order or larger in the case of 
the appearance probability, as discussed above. 

The question of how to interpret the total time t in the transition probabili- 
ties, along with the conceptual differences between appearance and disappearance 
(insofar as the oscillatory contributions), suggests a re-examination of how these 
probabilities should be calculated. 

It proves illuminating to understand the result above from "old-fashioned" time 
dependent perturbation theory with a finite time interval tf, fy. To this end, we will 
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need the second order matrix element 
A a ^f) = -g 2 J d 3 x 1 d 3 x 2 J dti J dt 



(W(k D )Jp(p D );l a (ps)\W(x u t 1 )l(xi,ti)v(xi,t 1 )W(x2 ) t 2 )l(x2,t2)v(x2^ 

(59) 

where we have suppressed the flavor indices in the interaction for simplicity of 
notation. We note that in the above expression time is ordered tf > t\ > t 2 . There 
are two Wick contractions corresponding to the processes displayed in Fig. ([3]). 
In diagram (J5^i) the charged lepton at the near detector Is is created at t 2 < ti, 
namely before the charged lepton at the far detector Id, whereas Fig. (J5}d) displays 
the opposite process in which the charged lepton at the far detector is created before 
the charged lepton at the near detector. 

= ~9 2 n (2tt) 3 5(k s - Ps -k D - p D ) e lAE * x 



U aj \ 2i sin (AM) e ~ iAE i - e -^t e i(E D +E S )i 



2%(?) | (E s - %(<z))A£ (n 3 (q) - B s )(fij(?) - E D ) 
4*La = ~9 2 n (2tt) 3 S(k s - p s -k D - Pd) <? ae * x 



U aJ j - 2i sin (AE±) 



e lAE i -e »>"e 



1^,(61) 

where q — ks — ps ■ It is straightforward to confirm that A^\a + A^\^ coincides 



2n j (q) \ (E s + Qj(q))AE (%(<?) + £fe)(fi,(g) + E D ) 

; — ps- It is straightforward to confirm 
with the results (|38l42j) at finite time taking S, e = 0. 



l a,S l (3,D l a,S l [3,D 




Fig. 3. The two contributions to the transition amplitude l|59|l from the two different Wick 
contractions. For both diagrams t± > ti- 



In process (a), the interaction vertex annihilates the Ws in the initial state and 
creates both a neutrino and the charged lepton lp g at the space-time coordin ates 
{x 2 ,t 2 ) in a superposition of product Fock states, namely in an entangled state^^i. 
In process (&), the initial Ws and the final lepton lp t s are annihilated and created 
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at respectively and the neutrino and Wd, 1<x,d are created at (a?2,£2), with 

the time ordering t\ > £2- The combination of both contributions yields the time 
ordered product as usual. However, for a long-baseline experiment when the size of 
the near and far detectors are much smaller than the distance between them, these 
processes should not contribute on equal footing: the charged lepton produced at 
the source (near detector) will be detected much earlier than the charged lepton 
produced at the far detector. 

A somewhat extreme example which can illustrate the point is SN1987a in the 
Large Magellanic Cloud. The charged leptons produced in the explosion along with 
the concomitant neutrinos are trapped, or "detected," in the optically thick medium, 
whereas the neutrinos produced at the source are detected at a distance ~ 51 kpc 
away from the production region. The detection of the charged lepton in the near de- 
tector disentangles the quantum state. This suggests that long-baseline experiments 
actually involve two times: the time at which the charged lepton produced with the 
neutrino is detected in the near detector and the time at which the charged lepton 
produced by the neutrino is detected in the far detector. Obviously, only the process 
described by the amplitude A^l+p above describes this physical situation. However, 
this amplitude does not describe the process of measurement of the charged lepton 
at the near detector, namely, the disentanglement of the charged lepton and the 
neutrino. The correct description requires addressing the issue of coherence and 
entanglement. 

4. Coherence, entanglement and oscillations 

In order to study aspects of coherence we consider a simplified interaction La- 
grangian density 

Ci = gW ev e — gW e(cos9 vi + sin 1/3), (62) 

focusing only on one lepton, which we refer to as the "electron" to simplify the dis- 
cussion. The full coupling as in Eq. can be treated similarly without modifying 
the main conclusions. Although W may be interpreted as a charged vector boson, 
the analysis is obviously the same if it describes a pion field. 

We can study aspects of coherence by focusing on the Fock state obtained upon 
evolution of the decaying initial state. We consider a plane-wave Fock initial state 
W(k)) at ti = 0. By time t, this initial state has evolved into |W(fe)) e~ tE £ 1 + 
^ e (i)). To lowest order in the interaction we find the second term to be 

\W(k)), 
(63) 

where all the fields are in the interaction picture. Though the field operator W can 
either annihilate the initial state or create another W particle, the state with two 



|*e(*)> =W< 



,//„/; / dl > 



d 3 3 



W(x, t') e(x, i')(cos 9 v x (x , t')+sm 6 u 2 (x, t')) 
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W particles features faster oscillations that will average out. In what follows, we 
consider only the Fock state resulting from the annihilation, leading to the state 



|*e(*)> 



2J2VEY 



9 



AW -Ei-a 3tlt )t 



1 



cos f 



E~ 



jifEf-Pf-flyJi _ j 

(Ef^Ei-n hP ) 



p = k — q 



(64) 



in which the electron and the neutrinos are entangled^ The neutrino state that is 
entangled with the muon is obtained from (|64[) by replacing cos — ¥ — sin 9 ; sin 6 — > 

COS0. 



4.1. Unobserved daughter particles: time evolution of the density 
matrix 



If the electrons (or daughter particle in Ref.(p2l)) are not observed, they can be 
traced out of the density matrix \^ e {t)) e (t)\. This gives the reduced density 



matrix 

Pr {t) = Tr e |# e (i))(# e (*)| 



E 



sm 



WE™ 

k q 



COS 



&2,p E~ 



Qi,i> E~ 
sin 26* 



Vl,p){vi,p 



V2,p)(v2,p 



E\ - fl2,p) I 1 1 2 



{Ef-m-^)/2 



sm i 



W 



(Ef -EZ-Sl ltff )/2 
S in((£f-£|-n 2 ,^) 1 [sin ((£ fc F - E| - n X) 

(sf -£|-n 2)# )/2 



V) 2 



+ e « 



Vl,p){v2,p 



p = k — q. (65) 



This expression contains remarkable information. The function sin 2 (a;t)/a; 2 is the 
usual "diffraction" function of Fermi's Golden rule. In the formal long time limit 
sin 2 (xt)/x 2 — > ntS(x), the first two terms of the density matrix, which are the 
diagonal entries in the mass basis, describe the production process of the mass 
eigenstates. As will be seen below, in the long time limit, the time derivative of 
these two terms yields the production rate for the mass eigenstates. In the formal 



b The result for the wavefunction in Ref . 1221 may be understood using a (non-perturbative) Wigner- 
Weisskopf approximation for the decaying parent particle, replacing Eyy — » Eyy — iFyy. Taking 
t 3> l/r^jr in the integral replaces the brackets in I I 6411 by l/(E\v — E e — fij —i^w) whose absolute 
value is proportional to S(E\y — E e — Qj)/Fyy. 



January 26, 2012 1:10 WSPC/INSTRUCTION FILE nuoscillation-IJMPA- 
revised 



20 Jun Wu, Jimmy A. Hutasoit, Daniel Boyanovsky, Richard Holman 

t — > oo limit, these are the diagonal terms obtained in Ref. ^321, where in that refer- 
ence, the product of delta functions is again understood as the total time elapsed 
times an energy conserving delta function. 



The off-diagonal terms in the last line of (|6"5)) describe the "coherences" and 
display the oscillatory phases from the interference of the mass eigenstates. The 
time dependent factors of the off-diagonal density matrix elements determine the 
coherence between the mass eigenstates and will be a ubiquitous contribution in 
the real time description of oscillations that follows below. The functions 



f±(x,t;A) 



2 sin 



• ± A) 



± A) 



Ef 



rip* 



A 



40 



(66) 



as functions of x are strongly peaked at x± A = with height t and width ~ 2n/t. In 
the infinite time limit f±(x, t, A) — > 27r <5(x±A) and thus, their product would vanish 
in this limit, leading to the vanishing of the coherence. This is the result obtained 
in Ref. ^2. However, at finite time t, they feature a non- vanishing overlap when 
2 A < 2%/t. We recognize this precisely as the condition for oscillations. We note 
that t ~ 7r/A yields a macro scopically large time scale. The functions f±(x,t, A) 
and their products are depicted in Figs. (|4I5[) for the values A = 0.1, t = 40, 100, 
respectively. 



< 



2 A 




; f+(x,40,0.1)/ \ / 


\f_(x,40,0.1) 







< 



_2A 




r l(x.lOO.O.I) 1 


lf_(x,100,0.1) 





,1 



-0.6 -OA 



Fig. 4. The functions f±(x, t, A) vs. x for t = 40, 100, A = 0.1 



It is straightforward to find 



/_(£,*, A)/+0M, A) 



sin(At) 



s'm[(x - A)t] sin[(x + A)t] 



A 

cos(At) 
A 



(x — A) 



(x + A) 



sm 2 [(x- A)|] sin 2 [(x + A)|] 



(x-A) 



(x + A) 



(67) 



In the long time limit, the terms in the first bracket yield a sum of delta functions at 
x = ±A, whereas the second terms are similar to a principal part. Upon integrating 
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Fig. 5. The products f-(x,t, A)f + (x,t, A) vs. x for t = 40, 100,A = 0.1 



the product of /-/+ with functions of compact support, the contribution from the 
second line in (|67|) is negligible in the long time limit. Therefore, the long time limit 
of /+/- can be replaced by 



/_(»,t,A)/ + (x,t,A)=7r 



sin(A t) 
A 



S(x - A) + S(x + A) 



(68) 



For a large time t <IC 2n/A, it follows that the product f-(x, i, A)f + (x, t, A) ~ 
7r i — A) + 5(x + A)] grows linearly in time but is bound in time. For t > 2ir/A, 
it oscillates with frequency 2n/A. Therefore, we conclude that upon integration 
with a smooth density of states, the off-diagonal terms in the density matrix grow 
linearly in time for t <C t osc = 27r/A, but feature a bound oscillatory behavior of 
frequency A for t >2tt/ A. 



Whereas the diagonal terms, i.e., the first two terms in the reduced density 
matrix (|65j) . arc proportional to 4sin 2 (xi)/x 2 — > 2n tS(x), the coherences or the off 
diagonal terms are linear in time and of the same order as the diagonal elements 
for t < t os = 2tt/ A, but are of 0(1/Af) and oscillate fast compared to the diagonal 
terms for t 3> t osc . This is similar to the phenomenon observed in the transition 
probability in the previous section. This behavior is displayed in Fig. (j6]), where as 
an example we consider a smooth density of states and the integral 

/oo 
e- x2 f+(x,t,A)f-(x 7 t,A)dx. (69) 
-oo 

The case A = describes either of the diagonal terms and is linearly secular in 
time. This figure clearly shows the slow oscillations for t > 1/A. 

Therefore, the approximation (|68j) is reliable in the long time limit and upon 
integration with functions of compact support. We see that for large t, but tA <C 1, 
/+ (x, t, A)/_ (x, t, A) -> [S(x - A) + S(x + A)] and for A ->■ 0, the product yields 
2irt8(x). 
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Fig. 6. The integral /(t, A) = /f^ e~ x f+(x, t, A)/_ (x, t, A)dx vs. t for A = 0, 0.1, 0.3 



The reduced density matrix (|65l) allows us to obtain the time evolution of the 
neutrino populations and coherences, namely 

n, (Jp, t) = Tr p r {t)a\ (p)a,i (p) ; C,j (p, i) = Tr p r (t)a\ (p)a,j (p) , i ^ j (70) 

where the annihilation and creation operators are in the Schroedinger picture. In 
the long time limit and using the results above we find 



m(p,t) = tT 1 (p) cos 2 6 ■ n 2 {p,t) = t r 2 (p)sin 2 6», (71) 



where 



are the partial widths for the decay of the W into the charged lepton and the 
neutrino mass eigenstates. Similarly, 



C 12 (p,t) = Cl(p,t) 



2irg 2 sm29 sinftA] _ iAt 



32SF A 

fe 



e 



(73) 



dm(p,t)/dt yields the production rate of the neutrino mass eigenstates from the 
dec ay of the W, and the coherences Cy are non- vanishing at any finite time. In Rcf. 
1221, these coherences vanish as a consequence of the product of delta functions on the 
different mass shells of the mass eigenstates. The coherences vanish in the formal 
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infinite time limit because of the oscillatory behavior averages out on time scales 
t 3> 1/A. However, we would like emphasize that on the experimental situation, 
the time scales involved (or rather length scales) are of order 1/ A as these are the 
scales on which oscillatory phenomena are revealed. Taking fl± ~ ~ H in the 
denominators in (172I73[) . it follows that 

„ , _ , sin2#sin[iAl „- At 1 r „ „ -, 

C 12 (p,t)~ — 1-L e -^-[T 1 + T 2 ] . (74) 

Therefore, the coherences are of the same order of the population terms on time 
scales t < 1/A, but average out for t ^> 1/A. This clearly shows that coherence is 
maintained over the oscillation time scale. 



4.2. Disentangling the neutrino: a two-time measurement 



As we have discussed previously, a long baseline experiment is actually a two time 
measurement, as the charged lepton produced at the interaction vertex at the source 
is detected by the near detector. This "measurement" of the charged lepton disen- 
tangles the neutrinos from the charged lepton in the quantum state (|64|) ^21. The 
detection of the charged lepton at the near detector projects the quantum state (|64|) 
at the observation time ts onto the single particle charged lepton state e 
disentangling the neutrino states into the "collapsed" state 

\V e (P,t s )) = (e d \*(t s )) e lE * ts 



-iEtts 



e o) 




E s = Ejy -E%; P = k-Q 
We note that up to a phase the coefficient functions that multiply \fi,2) are 



(Es 



(76) 



respectively. In the limit of ts —> oo, these become 2tt8(Es — Sli) and 2ir5(Es — O2), 
respectively. Therefore, in this limit, for a fixed Es, one of the quantum states will 
be projected out. However, as we insist on keeping a finite time interval, we will 
keep ts finite. 

The state V e (P,ts)) then evolves forward in time with the full Hamiltonian 

\V e (P,t)) =e- iHot U(t 7 ts)e tHots \V e (P,t s )). (77) 
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The "free" evolution is obtained by setting to lowest order U(t,ts) — 1, leading to 



\V e (P,t))=i 




The phase factors e~ in i ts ' 2 multiplying each mass eigenstate are the consequence 
of the phase build-up during the time evolution from the production vertex until 
the detection of the charged lepton and the collapse of the wave function. These 
can be absorbed into the definition of the states \i/x 2} at the fixed time ts- 



(78) 



The expression (|78|) features the factors 



{Es-^/2 ' (79) 

which as ts — » 00 becomes 2ttS(Es — These factors, which are a direct conse- 
quence of the neutrino state being produced by the decay of the "parent" particle 
(here the W) into an entangled quantum state, distinguish Eq. (|T8[) from the famil- 
iar quantum mechanical description. These factors emerge from the (approximate) 
energy conservation at the decay vertex. Again, in the ts — > 00 limit, if the energy- 
momentum of the parent particle and the charged lepton are both certain, only 
one of the mass eigenstates will be produced but not both. However, writing VL12 
as in Eq. (|49]). it follows that for t s < -Ifi/Jm 2 , the width of the "diffraction" 
functions is much larger than the frequency difference A and there is a substantial 
overlap between these "approximate" delta functions. Only for ts > t osc = 2ir/A 
are the two peaks at Es — O = ±A actually resolved, whereas for t <C t osc , the two 
peaks are unresolved, "blurred" into one broad peak at fl Thus, we can use the 
approximation 



sm 



(E s - fii) 



(Es 



(Es 



(E s 

for t s < 2tt/A. 



fii)/2 



{E S - Sh)/2 



(E s - n)/2 



(80) 



To illustrate the validity of the above approximation, let us consider the case 
in which the typical size of the near detector is a few meters across. In a typical 
experiment, the charged lepton emerging from the interaction vertex travels this 
distance within a time scale ts ~ 10 -8 s, leading to an energy uncertainty above 
AE - h/t s - 10" 7 eV. Taking as an example 5m 2 - 10~ 4 eV 2 ; IT - E s - 100 MeV, 
it follows that Sm 2 /il ~ 10 -12 eV AE. Therefore, the detection at the near de- 
tector cannot resolve the energy difference between the mass eigenstates and the 
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approximation (1801) is justified. 

Another approximation we can use in (1751) is Q,\ ~ — Ct, for SI ^> A. Ab- 
sorbing the phase factors e~ ir2 j's/ 2 m t definition of the states the time 
evolved disentangled state is then approximately given by 

sine \v 2 ? ) e- m ^ {t - ts) +cos6 \v l P ) e-^i.f 

(81) 

for t s < 2tt/A. 

The state inside the brackets is identified as the usual quantum mechanical state 
that is time evolved from the "electron" neutrino state, which is prepared initially 
at ts- From this analysis, we see that there are two conditions required for the 
disentangled neutrino state to be identified with the familiar quantum mechanical 
state. The two conditions are Sm 2 /CI -C 1 and ts <C t osc ~ 2irQ,/8m 2 . The former 
is always satisfied for neutrinos with 5m 2 ~ 10~ 3 — 10~ 4 eV 2 , il > few MeV, 
while the latter is fulfilled for near-detection of the charged lepton in long-baseline 
experiments. The latter condition implies that the neutrino state is disentangled 
before oscillations can occur. In a long-baseline experiment this is achieved if the 
charged lepton, which is entangled with the neutrino at the production vertex, is 
measured at the near detector . 



\Ve(P,t)) 



>fj 



4.3. Transition amplitudes and event rates 



The number of charged lepton events with momentum Q at the near detector, at 
time ts is given by 

n e (Q,t s ) = (^e(t s )\at(Q)a e (Q)\^ e (t s )) = (V e (P,t s )\V e (P,t s )). (82) 

For ts -C t osc = 2n/A and ft ^> A, using the approximations leading to (|8T|) . we 
obtain the differential detection rate at the near detector 



(2tt) s 



dN, 



(ND) 



d 4 xd 3 Q 



dn e (Q,t s ) 
dt s 

2.9 2 

8veWe% n 

k Q 

2ng 2 _ 

sveWe^ n 



sm 



(E s - Sl)t s 



{E s - (I) 

s(E s -n), 



(83) 



where at the last step we have replaced the diffraction function by the delta func- 
tion. This can be justified as follows. For ts ~ 10 -8 s, the width of this function 
(the resolution) in energy is ~ 10" 7 eV. Since the typical energy in a long-baseline 
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experiment is > 40 — 100 MeV, the error incurred in replacing the diffraction func- 
tion by a delta function is smaller than one part in 10 15 . 

We can also obtain the transition amplitude for the disentangled state to produce 
a final charged lepton and another W particle at the far detector at time t d , where 
try — ts ~ L and L is the baseline. It is given by 

A a ^ p = (W(kD),h(pD)\e- iH ^- ts ^ \V e {P,t s )) 

= e - iEDtD {W(k D ),h(p D )\ U(t D ,t s )e iHots \V e (P,ts)}- 

The disappearance and appearance amplitudes are then given by 

Ae-> e = -g 2 n (2%) 3 S(k s -ps~k D - Pd) x 



(84) 
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'(Es-n 2iP )± 



sm 
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(85) 



and 



-g 2 U(2nf5(ks - ps ~ k D - p D ) 



-iCt. s- 
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i£s — f^i 5 



i,p/ 2 



sin26> 
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[E D -Sl lP ) (t D -t s )/2 
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(E S - n 2 «)/2 



(£ D - Q l i3 )/2 
(E D -n 2 p)(t D -ts)/2 



(E D - n 2 p)/2 



(86) 



with P = kg—pg. In these expressions, we have used the same notation as in section 
(j3|), where LT is given by (|22j) . Implementing the same approximations leading to 
the factorized state (j5Tj). namely Q, ^> A and ts A <C 1, we find the disappearance 
and appearance probabilities 

V e ^e(t D ) = (2ir) 3 V 5(k S - PS - k D - PD)2TTt S 5(E S - Tip) 

I cos 4 6 fl (x, t, A) + sin 4 6 f 2 (x, t, A) + 2 cos 2 sin 2 6 cos(tA)f + (x, t, A)/_ (x, t, A) 



P. 



e—tfj. 



(t D ) = (§^) V) 3 V 5(k s - pg -k D - p D ) 2n t s S(E S -Up) 



fl(x, t, A) + fl (x, t, A) - 2 cos(t A)/ + (a, t, A)/_ (x, t,A) >, 



(87) 
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where t = to — ts and x = Ed — Qp- Here 
and f± are given by Eq. 1661 



E fs , E D -E % 



E PD> 



In the long time limit, using f± (x, t, A) — > 2ir t S(x ± A) and we find 



Pe^e(t D ) = (^)\2n) 5 VS(k s -p s -k D -p D ) t s 5(E s ~Qp) 



cos 4 9 t 6(x + A) + sin 4 9 t 6(x - A) + 2 cos 2 9 sin 2 9 sm (^ A ) i 

sin 2 2(9 



<5(x + A) +<5(a; - A) 



V e ^u(t D ) = (^-) 2 {2nfV8{k s -p s -k D ^p D ) t s S{E s -n p ) 



tS(x + A) +tS(x - A) 



sin(2iA) 1 
2A 2 



5(x + A) + 5(x - A) 



These transition probabilities are very different from the ones obtained in sec- 
tion p.2j) and from those obtained from the S-matrix approach. They feature the 
two time scales t$ and t = tjj — ts associated with the measurements at the near 
and far detector. They also feature energy conserving delta functions associated 
with the different mass eigenstates. 



(89) 



There is a further simplification when Q 3> A. In this regime, when the probabil- 
ities (|89I89[) are integrated over a smooth density of states, the delta functions corre- 
sponding to the mass eigenstates yield the density of states at values Ed = ft T A. 
In typical experiments, where f2 ~ 100 MeV and Sm 2 ~ 10~ 3 eV 2 , the density of 
final states must vary dramatically near f2 to resolve the small interval A, with 
A/fl < 10 -19 . Therefore, understanding the probabilities as being integrated with 
a smooth final density of states insensitive to the mass difference, we can approx- 
imate 5(x ± A) ~ S(x). In this case, we can approximate the above expressions 

by 

V e ^e{t D ) = (ff^;) 2 W) 5 V S(k s - ps - k D - pd) t s 5(Es-Up)5(E D -Tip) 



t 



5 2 ri\ 2 



2 cos 2 9 sin 2 9 
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2Q p 



T e ^u(t D ) = I ~= ) {^T V S(k s - ps - k D - p D ) t s 6(E s -np)S(E D -np) 



(90) 

^ , sin 2 29 



t 



sin(2iA) 1 



2A 



(91) 



The product S(Es — Qp) S(Ed — Qp) is an approximate energy conservation at both 
production and detection vertices, where we have neglected A, which is twice the 
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energy difference between the mass eigenstates. However, even under these justified 
approximations, the probabilities (|90I91I) are very different from those obtained by 
the S-matrix approach, even after including the finite time corrections discussed 
in section (13.2p . They also differ greatly from the transition probabilities from the 
simple quantum mechanical argument. 

Further insight can be gained by obtaining the phase space distribution of the 
number of charged leptons I = e, p at the far detector 

(2^) 3 =MPD,t D ) = (V e ($,t D ,ts)\<$(pD)ai(p D )\V e ($,t D ,ts)). (92) 

Here, 

\V e (Q,t D ,t s )) =e- zH ° tD U(t D ,t s )e lHots \V e (Q,t s )) (93) 

is the neutrino state disentangled at ts at the near detector and time-evolved until 
it is detected at the far detector at time to- Not surprisingly, since up to order g 2 , 
the time evolved state contains a single lepton Fock state, we find that 

(2tt) 3 -^§^=7W^) ; {2*Y-^^=V e ^(t D ), (94) 
arxarpjj a^xd^pu 

with the probabilities P e ->e(*z>) and V e -yn(tD) are given by (|9QI9ip . 

Taking the time derivative with respect to to, we obtain the differential charged 
leptons event rates at the far detector 

{27r)3 d^xdtd^D = ( 2 *f v5 $s-Ps-k D -p D ) t s S(E s -Q p )S(E D -Qp) 

\ cos 4 9 + sin 4 6> + 2 cos 2 6 sin 2 6 cos(2i A) \ , (95) 



3 ^x%r D = (& {27T)5v6Cks ~ " s - %D - pD) ts s{Es - n ^ s{Ed ~ ^ ^ 



jl -cos(2tA)|. (96) 

Remarkably, these rates can be simply written as 

, dNi D , dN ND 

^d^Jd^ = ^ 3 d^ P ^ dT ^> < 97 > 



where we have used the expression (|83|) for the differential charged lepton event rate 
at the near detector, integrated in ts, 

&T V ^ WP = ^y^/j^ S $s - Ps -k D - p D ) S(E d - U p ) (98) 

k D PD 
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is the charged lepton production rate from the reaction v$ — ¥ Wlp for a flavor 
neutrino of energy f2 and P a ->p(t) are the disappearance (a = fJ) or appearance 
(a 7^ j3) quantum mechanical transition probabilities. 

The remarkable aspect of the final result (|98[) is the factorization of the different 
processes contributing to the far detector event rate, namely the number of events 
at the near detector multiplies the quantum mechanical transition probability which 
in turn multiplies the production rate at the vertex in the far detector. This factor- 
ization is a distinct consequence of the two time analysis, of the disentanglement 
of the neutrino at the near detector, along with the approximations invoked in the 
resolution of the energy conserving delta functions. We emphasize that this factor- 
ization in terms of the quantum mechanical transition probabilities only applies to 
the detection rate defined by taking the time derivative, but not to the total number 
of events or to the rate defined by simply dividing by the time scale. 



5. A model for the GSI anomaly 

Recent experiments at the Experimental Storage Ring (ESR) at GSI in Darm- 
stadt have revealed an unexpected time dependent modulation in the population of 
daugh ter ions 140 Ce 58+ from the electron capture decay 140 Pt- 58 + — > 140 (7e 58+ + v e 
| 68 | 69 |^ a phenomenon that has been dubbed the "GSI anomaly." There have been 
some works that try to explain this remarkable time dependent rate of the parent 
ion dec ay as an interference be tween the neutrino mass eigenstates in the decay 
reaction | 68 | 69 | 70 | 71 | 72 | 73 | 74 | 31 [ However, this interpretation has been re-examined 
and criticized E3E 4 | 65 | 66 | 67 | 

The authors in Refs. | 70 | 71 | 72 | 73 | 74 | obtain the decay rate of the parent ion by 
adding coherently the amplitudes and then obtaining the probabilities, in which case 
the modulation arises from the interfere nce of the m ass eigenstates in the squaring 
of the amplitudes. The authors in Refs. I 63 | 64 | 65 |66| cr jticize this approach, stating 
that it is not the amplitudes that must be summed coherently but the probabilities, 
corresponding to an incoherent addition of the contributions from the mass eigen- 
states. This approach does not lead to any modulation as the amplitudes for the 
different mass eigenstates do not interfere. 

Rather than focusing on either one of these approaches, we analyze the situation 
differently, by obtaining the time evolution of the population of the parent 
and daughter particles. In this approach, we recognize directly the decay rate of 
the parent particles (or production rate of the daughter) without the necessity to 
invoke a coherent sum over amplitudes or a sum over probabilities. 

We study the time evolution of the populations by modeling the situation in 
terms of the decay of a parent particle via charged current interactions into a charged 
lepton and its associated neutrino. The interaction vertex is given by W ev e , where 
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we simply take the parent to be the W and the daughter to be the e-charged lepton. 



Let us consider an initial parent particle state |W(A;)) that is prepared at time 
t = 0. The evolution of the number of parent and daughter particles is obtained 
from 



N w {k,t) = (W(k)\ e tHt ok(*W(£) e~ im \W(k)), 
n e (Q,t) = (W(k)\ e lHt at(Q)a e (Q) e~ lHt \W(k)), 



(99) 



where the annihilation and creation operators are in the Schroedinger picture. We 
note that e~ lHt = e~ lHat U(t,0) and that the number operators commute with the 
free field Hamiltonian. We also note that U(t, 0) | W) = \ W) + | * e (t)) (1) + 1 * e (i)) (2) + 
• • • where 



y e (t)) (1) =ig I dh 
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(100) 



and 



|*e(*)) 



(2) 



dh / d 3 X! 
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W(xi, h)e(xi, ti)u e (x i, ti)) 



W(x 2 ,t 2 )e(x 2 ,t 2 )v e (x 2 ,t 2 )) |W(fc))Gl01) 



with v e = cos#i/i + sin#zy 2 . Since the (entangled) state \^ e (t))^ has one daughter 
particle (electron) and the initial state has none, it is clear that to lowest order, the 
number of daughter particles is 



n e (Q,t) = W(* e (i)|a|(Q)a e (Q)|*e(t))« = (V e (Q,t)\V e {Q,t)), 



(102) 



where \V e (Q,t)) is the "collapsed" state (|75]) . We find the production rate of the 
daughter particle 



dn e (Q,t) 
dt 



8V E^y El 

k Q 
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If at this stage we take the long time limit, using the identity 
2 sin [(if -E^-U 3 )t] 



2tt5(E? -E^-Qj), 



(eW - E e - - n 2 

y k Q 

(103) 



(104) 



we find the total number of daughter particles produced as a function of time 
n e(t) = J2q n e (Q, t) is given by 



n e (t) 



Li cos 2 9 + T 2 sin 2 i 



(105) 



where Ti^ are the partial widths given by (|72p . In the rate (I103p . there is no interfer- 
ence between the mass eigenstates since |fi, 2 ) are orthogonal. It is straightforward 



January 26, 2012 1:10 WSPC/INSTRUCTION FILE nuoscillation-IJMPA- 
revised 



Neutrino Oscillations, Entanglement and Coherence 31 

to confirm that 

n e (t) = W<¥ e (t)|iMt)> (1) = ^2(Ve(Q,t)\vMt)) » (106) 

Q 

a result that will play an important role below. 

The calculation of the parent population is slightly more involved. Since the 
first order state does not have any parent particle W, we must consider the second 
order state (|101[) . To second order, there are several contributions, but the only one 
that is relevant is the process in which the first vertex at (x2,t2) annihilates the 
initial W creating the intermediate state with one (e, u e ) entangled pair, while the 
second interaction vertex at {x\,t\) annihilates this (e,v) pair in the intermediate 
state and creates the W, which has non- vanishing overlap with \ W). This process 
is depicted in Fig. ([7]) and is recognized as the self-energy of the parent particle. 




Thus to lowest order in <?, 

N w {k,t) = l + 2Re[(VF|* e (t)) (2) . 

It proves more convenient to calculate dNw(k,t)/dt, for which we 

dN w (k,t) _ 2g 2 f d 3 Q 
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(108) 



c Effcctively, we are obtaining the Boltzmann equation for the parent particle, neglecting the build- 
up of the population. 
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In the long time limit this becomes 

dN w {k 7 t) _ 
dt 



I', cos 2 9 + T 2 sin 2 i 



(109) 



Clearly, dNyy(k, t)/dt = —dn e (t)/dt as the decay of the parent population results in 
the growth of the daughter population with the same rate. This is a consequence of 
unitarity and we can see this by substituting (I106I107|) into the unitarity condition 

(W(k)\u\t,0)U(t,0)\W(k)) = 1 = l+ (1) (* e (t)|* e (t)) (1) +2Re[(t¥|* e (i)) (2) )j+C(5 3 )- 

(110) 

The interpretation of (|109j) in terms of the Feynman diagram ([7]) also makes the 
conclusion of lack of interference manifest. The decay rate is the imaginary part of 
the self-energy. Since the correct propagating degrees of freedom are the neutrino 
mass eigenstates, the self-energy is the sum of the self-energy diagrams with the 
neutrino mass eigenstates propagating inside the loop. Therefore, it simply follows 
that the total decay width is the sum of the partial decay widths on the mass eigen- 
states without interference. The real time calculation presented above confirms this 
result directly from the evoluti on of the parent and d aughter populations. Thus, 
we confirm the analysis of Refs. | 41 | 4z | 4j | 44 | 4o | 4o | oo | o l \ there is no interference 
of mass eigenstates and we conclude that the GSI anomaly cannot be explained in 
terms of the interference of mass eigenstates in the decay. 



6. Conclusions and Discussions 

In this article, we carry out an in-depth review on our recent discovery on the the- 
ory of neutrino oscillations, we study the dynamics of mixing and oscillations in 
quantum field theory directly in real and finite time. The setting is a bosonic model 
that reliably describes charged current-weak interactions. This allows us to extract 
the relevant aspects without the peripheral complications associated with spinors. 
We begin by obtaining the in-out transition amplitudes and probabilities in a long 
but finite time interval for appearance and disappearance processes and compare 
these to the S-matrix results. To illustrate the effects of the finite time interval, we 
consider the simplest setting with plane wave states, discussing how a wave-packet 
treatment needs to be modified to include the finite time contributions. We find 
finite time contributions that display the oscillatory behavior resulting from the 
interference of mass eigenstates and show that these corrections can be of the same 
order of or larger than the S-matrix result for appearance experiments. A deeper 
analysis of the different contributions lead us to argue that the in-out treatment 
is ill-suited to describe long-baseline experiments. The main reason is that in the 
production vertex at the source or near detector, the neutrino is produced in an 
entangled state with the charged lepton. 



January 26, 2012 1:10 WSPC/INSTRUCTION FILE nuoscillation-IJMPA- 
revised 



Neutrino Oscillations, Entanglement and Coherence 33 

The concept of the hnite time interval is compatible with the recent observation 
that the neut rino forms a n entangled state with the charged lepton produced at the 
near detector | 46 | 32 | 56 [57| anc j ena bi e us to study the details of the dynamics of this 
entanglement. If the charged lepton (or daughter particle) produced at the near 
detector is not measured, tracing out this degree of freedom yields a density matrix 
for the neutrino. The off-diagonal density matrix elements in the mass basis are a 
measure of coherence and show distinctly the oscillatory interference effects. Diag- 
onal and off diagonal matrix elements are of the same order of magnitude during 
time scales t < t osc ~ 2E/Sm 2 , with E the typical neutrino energy, revealing that 
coherence survives during the time scale t OS c- On the other hand, the measurement 
of the charged lepton at the near detector disentangles the neutrino state, and it 
is the further time evolution of this disentangled state with the interaction Hamil- 
tonian that leads to the production of charged leptons at the far detector. Thus, 
the process of production and detection in long-baseline experiments involves two 
different time scales: the measurement of the charged lepton at the near detector 
determines the first time scale at which the neutrino state is disentangled, while 
the measurement of the charged lepton at the far detector is the second and longer 
time scale. 

The form of the amplitudes (|53I54|) and probabilities (|55|56|) obtained from the 
in-out formulation are fundamentally different from those obtained from the time 
evolution of the disentangled state ([551 I56"[) . We can see why this is so by noticing 
the difference in how and when the interaction vertices act on the state. In the 
transition amplitudes ([531541) . the second order matrix element involves the time 
ordered product or alternatively the nested time integrals (see ) , which result in 
overall energy conservation. On the other hand, in the amplitudes obtained from the 
disentangled state, the perturbation acts at different and non-overlapping times. A 
first order vertex creates and evolves the entangled state until the measurement of 
the charged lepton disentangles the neutrino state, while another first order vertex 
propagates the disentangled neutrino and creates the final charged lepton that is 
measured at the far detector. There is another way to understand this. As made 
explicit by time dependent perturbation theory (see subsection 13. 3[) . there are two 
contributions (corresponding to the time ordering) whose sum yields the in-out am- 
plitudes. These contributions are also manifest in the in-out probabilities. Whereas 
in the second case, there is only one contribution (as discussed in subsection 13. 3[) . 
These differences explain the fundamentally different form between the amplitudes 
and probabilities obtained from the in-out formulation, even at finite time, and 
those obtained from the disentanglement at the near detector and further evolution 
of the neutrino that leads to the charged lepton at the far detector. 

For long baseline neutrino oscillation experiments corresponding to the case of 
the disentangling time much smaller than the baseline, We obtain the charged lepton 
event rate at the far detector and factorize this rate into a product of the number 
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of charged lepton events per phase space at the near detector times the quantum 
mechanical transition probability, times the differential production rate at the far 
detector. This factorization is a direct consequence of the disentanglement of the 
neutrino state at the near detector along with the approximation that the density 
of final states is a smooth function and for large energies, it is not sensitive to the 
energy difference of the mass eigenstates. This factorization was assumed in Ref. ^1 
and physical arguments were proposed for its validity in Ref. ^1 However, in this 
article, we show that the time evolution of the disentangled neutrino state unam- 
biguously leads to the factorization, under the approximations discussed above. We 
would like to emphasize that the factorization in terms of the quantum mechanical 
probabilities is only valid for the event rate at the far detector. It is not valid for the 
total number of events. This latter quantity is also factorizable, but not in terms of 
the transition probability, but its time integrated version, which leads to a different 
energy dependence of the oscillatory terms as can be gleaned from the second lines 
in Eqs. (j90l9T|) . 

More interestingly, the appearance and disappearance amplitudes (|85l l86|) de- 
pend on both the disentangling time scale and the final detection ti me scale. For 
short baseline neutrino oscillation experiments, e.g., LSND experiment I 81 | 82 | 83 | 84 |^ 
the disentangling time can be of the same order as the detection time, which is 
also of the same order as the oscillation time scale. In this case, the disentangling 
time introduces an extra modulation with energy (see (|851 I86p ) which may yield 
phenomcnologically interesting modifications in the interpretation and analysis of 
data. Based on the dynamics of en tanglement and disentanglement discussed here, 
a recent study by Boyanovsky analyzes the influence of the disentanglement on 
neutrino oscillation formulas of short baseline neutrino oscillation experiments by 
introducing yet another time scale describing the decay length of the parent particle. 
It was found that both the disentanglement length scale and the source lifetime lead 
to a suppression of the oscillation probabilities in short baseline experiments, and 
imply that fits to the experimental data based on the standard quantum mechani- 
cal formulas underestimate both the mixing angle and the mass difference between 
different generations of neutrinos. 

Although our study has been carried out in terms of plane waves, we can ex- 
trapolate our results to include the case of wave packets. Let us consider the initial 
state in (|63|) . Instead of being the plane wave Fock state |VF(fc)), let it be the lo- 
calized state \W(S)} given by (fTTjl. Let us also assume that the charged lepton 
measured at the near detector is measured with a wave-function represented by a 
wave-packet (fTTj) \l e (S)). The neutrino is therefore disentangled into a wave-packet 
state, obtained from (|75[) by the convolution with the wave-packets of the W and 
e. This resulting state will be described by a wave packet whose localization length 
is determined by the production and disentangling processes. Upon further time 
evolution, this wave packet will spread and split into the two mass eigenstates and 
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produce the final charged lepton and W. The two mass eigenstates are also going to 
be described as wave packets. It is only when these wave packets begin to overlap 
with the far detector that the charged lepton at the far detector will be measured. 
The corresponding amplitudes will be non-vanishing when tp — ts is of the order 
of the baseline, with uncertainties of the order of the size of the far detector. These 
mildly localized wave packets lead to both momentum and energy uncertainty of 
the same order AE ~ Ap ~ h/a ~ K/ts ~ 10~ 7 eV. These uncertainties in energy 
and momentum are much larger than typical values for the neutrino energy dif- 
ferences 5m /fl ~ 10~ 12 eV. The issue of entanglement in neutrino oscillations was 
discussed within the wave packet approach in the framework of quantum field theory 
in refpH, which argued that S-matrix calculations give rise to the correct formulas 
of neutrino oscillations under realistic experimental setups. Indeed, the localization 
of wave packets leads to results somewhat equivalent to finite time calculations, and 
may yield standard formulas of neutrino oscillations for most realistic experimental 
setups. Both the spatial and time localization of the wave packets introduced in this 
reference are tantamount to a finite time description, although, as mentioned above 
the intermediate stage invoking the S-matrix actually refers to an in-out matrix 
element where the initial and final time are taken to -co, oo respectively. However, 
the interplay between the disentangling time and the final detection time and their 
influence on neutrino oscillation formulas are not directly captured by a wave packet 
analysis and, as we pointed out are fundamentally important and merit a thorough 
understanding. 

The results discussed above led us to study the GSI anomaly by directly ob- 
taining the decay rate of a parent and production rate of a daughter particle, by- 
passing the issue of whether amplitudes or probabilities for mass eigenstates must 
be summed. We show that these rates do not feature oscillations arising from the 
interference of mass eigenstates in the final state. We provide an alternative field 
theoretical explanation in terms of the imaginary part of the self-energy diagram of 
the parent particle. 

In this article, we have introduced the two time measurement approach by con- 
sidering the simplest model for neutrino oscillation. It would be interesting to apply 
this method to the mo re sophisticated models. A particularly interesting model is 
the 2 + 1 model of which features an unparticle sterile neutrino, the unsterile 
neutrino, along with the two active ones. The unparticle nature of the unsterile 
neutrino results in four momentum dependent mixing angles and non-trivial spec- 
tral densities of the neutrino mass eigenstates. Therefore, the quantum mechanical 
approach to the model of ^1 can on ly b e thought of as a proxy description and to 
understand its oscillation dynamics, the quantum field theoretical two time mea- 
surement approach is necessary. 
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